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Abstract 

Recently Beaton, de Gier and Guttmann proved a conjecture of Batchelor and Yung that 
the critical fugacity of self-avoiding walks interacting with (alternate) sites on the surface of the 
honeycomb lattice is 1 + \/2. A key identity used in that proof depends on the existence of a 
parafermionic observable for self-avoiding walks interacting with a surface on the honeycomb 
lattice. Despite the absence of a corresponding observable for SAW on the square and triangular 
lattices, we show that in the limit of large lattices, some of the consequences observed for the 
honeycomb lattice persist irrespective of lattice. This permits the accurate estimation of the 
critical fugacity for the corresponding problem for the square and triangular lattices. We consider 
both edge and site weighting, and results of unprecedented precision are achieved. We also prove 
the corresponding result for the edge-weighted case for the honeycomb lattice. 

1 Introduction 

Self-avoiding walks (SAW) in a half-space, originating at a site on the surface, are well-known and 
useful models of polymer adsorption, see [HI [21] for reviews. It is known [20] that the connective 
constant for such walks is the same as for the bulk case. To model surface adsorption, one adds 
a fugacity y = e" to sites or edges in the surface. Let c^(m) be the number of half-space walks 
of n-steps, with m monomers in the surface, and define the partition function as 



Zn{a) = E c:(m)e" 



m=0 



with a = -e/ksT, where e is the energy associated with a surface site (or edge), T is the absolute 
temperature and fcs is Boltzmann's constant. If e is sufficiently large and negative, the polymer 
adsorbs onto the surface, while if e is positive, the walk is repelled by the surface. It has been 
shown by Hammersley, Torrie and Whittington 10 in the case of the c?-dimensional hypercubic 
lattice that the limit 

lim n"^logZ„(a) = K{a) 

exists, where K{a) is the reduced, intensive, free-energy of the system. It is a convex, non- 
decreasing function of a, and therefore continuous and almost everywhere differentiable. Their 
discussion and proof apply, mutatis mutandis to the honeycomb and triangular lattices. 
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Figure 1: The figure sliows the two types of surface sites on the honeycomb lattice as indicated by 
solid and empty circles. A SAW starting at a and finishing at z is also shown. 



In the case of the honeycomb lattice there are two types of surface sites, marked as solid and 
empty circles in Figure [TJ hr most studies just one of the the two types is weighted; namely, 
those marked with solid circles in Figure [1] In this study we also allow for a surface weight on 
the second type of sites and we study the case where all surface sites carry the same weight. 

For a < 0, K{a) = log/i i20j, where fi is the connective constant for SAW on the given lattice. 
For a > 0, 

K{a) > max[log^, a]. 

This behaviour implies the existence of a critical value ac, such that, for the hyper-cubic lat- 
tice, < ac < log/i. The situation as a oo has only recently been rigorously established by 
Rychlewski and Whittington [19) . who proved that k(q;) is asymptotic to a in this regime. As 
illustrated in Figure (TJ it is convenient to attach weights y to only half of the sites along the sur- 
face to allow for simplihcations later on. In this case the bounds on ac become < ac < 2 log/i, 
or equivalently 1 < yc = e"" < fi^ . 

Various other quantities exhibit singular behaviour at yc- For example, the mean density of 
sites in the surface is given by 

^ ^ 1 Z,n"^<("^)/" ^ 1 dlogZ„{a) 
n ZmCi('^)y™ 9a 

In the limit of infinitely long walks one has 

dnja) 

^^"^ = 

From the behaviour of k given above, it can be seen that p{a) = for y < yc and p{a) > for 

y>yc- 



2 An identity for the honeycomb lattice with a boundary 

In a recent paper Beaton, de Gier and Guttmann j2j generalised a finite lattice identity by 
Duminil-Copin and Smirnov for the honeycomb lattice to the case where weights are intro- 
duced on alternating sites along a boundary as represented by the solid circles in Figure [TJ This 
resulted in a proof of a conjecture of Batchelor and Yung [1] that the critical surface fugacity 
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of self-avoiding walks interacting with (alternate) sites on the surface of the honeycomb lattice 
is 1 + ^/2. Here we briefly outline the results. 

Let H be the set of mid-edges on a half-plane of the honeycomb lattice. We define a domain 
i7 c iJ to be a simply connected collection of mid-edges. The set of sites adjacent to the mid- 
edges of is denoted V{^). Those mid-edges of VI which are adjacent to only one site in V{Vl) 
form 9f2. Since surface interactions are the focus of this article, we will insist that at least one 
site of V{Vl) lies on the boundary of the half-plane. 

Let 7 be a self-avoiding walk in a domain Q,. We denote by €(7) the number of sites occupied 
by 7, and by 1^(7) the number of contacts with the boundary. Define the following observable: 
for a 6 917, z e il, set 

F(a,z;x,?/,a):=F(z)= ^ Q--^<'W(^(a^z))^i(^)y.(i) ^ 

7(a^j2:)cf2 

where the sum is over all configurations 7 c for which the SAW goes from the mid-edge a to 
a mid-edge z. 1^(7(0 z)) is the winding angle of that self-avoiding walk. See Figure [T] for an 
example - the SAW shown there starts on the central mid-edge of the left boundary (shown as 
a) and ends at a mid-edge z. As the SAW runs from mid-edge to mid-edge, it acquires a weight 
X for each step and a weight y for each contact (shown as a solid circle) with the right hand side 
boundary. 



T 




Figure 2: Finite patch of the honeycomb lattice with a boundary. 



We define the following generating functions: 

St,l(x,2/):= E x'^'^^v''^^\ 

where the sums are over all configurations that have a contour from a to the a, /? or e,e 
boundaries respectively. For the SAW model in the dilute regime, the result proved in 2J for 
the n-vector model simplifies (in the case n = 0) to: 
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— j Ar.LixcV) + cosy- jET,L{xc,y) + _ ^^ Br^Lixc^v)- (1) 

where 

A simple corollary of (HJ is that at y = we have 
Corollary 1. 

1 = cos ^ j At,l (a;c , y* )+ cos ^ i;T,L (a^c , 2/* ) • (2) 

The importance of this result is that the generating function Bt,l for y = y* has disappeared 
from In [2] 'we proved that the critical surface fugacity j/c is equal to y* . Using this result 
and taking the limit L ^ oo, the geometry becomes a strip of width T, and the corollary then 
becomes ^ 

1 = cos^-^j AT(xc,yc) +cos^^ji;T(a;c,yc)- (3) 
In [1], we also proved that Et{xc, y) = for all < y < yc- So ^ simplifies further to 

1 = cos 



(^)AT(xc,yc). (4) 



This is a remarkable equation in that it implies that yc can be identified from the generating 
function AT{xc,y) for any width T, simply by solving equation Q. To show the power of this 
observation, note that virtually by inspection one can write down the solution for strip width 
0, which is 

Ao{x,y) = -. 

1 - a;^ 



Solving 1/ cos(^) = ylo(a;c,y), recalling Xc = l/\/2 + \/2, gives y = yc = 1 + \/2, the exact value 
of the critical fugacity. 

For other lattices, and even for the honeycomb lattice with interactions at every surface site, 
we do not have an equivalent identity, such as 1 = CaAT{xc,yc)- However if one plots AT{xc,y) 
versus y in these cases, one might be forgiven for thinking that such an identity exists. In Figure 
[3] we show a plot of ^7^(xc,y) versus y, for a range of strip widths T. To graphical accuracy it 
appears that there is a unique point of intersection for plots corresponding to higher values of 
T. Even finer resolution, see inset, suggests that this is the case. The actual small deviation 
can be seen from the data given in Table ID 

We denote by yc{T) the point of intersection of AT{xc,y) and ^T+i(2;c,y)- We observe that 
the sequence {yc{T)} is a monotone function of T. We argue, as in |3i, that in the scaling limit 
all two dimensional SAW models are given by the same conformal field theory. Since it is known 
for one of these models (i.e. honeycomb lattice SAW with alternate site interactions) that the 
critical point can be found by requiring certain contour integrals to vanish, it follows that in the 
scaling limit the same should be true for all two dimensional This is entirely consistent 

with our observations, and implies that limy^oo yc{T) = yc- 

This then suggests a potentially powerful new numerical approach to estimating y^. One 
calculates the generating functions Axixcy), for all strip widths T = 0, 1, 2, . . . Tmaxj uses these 
to calculate yc{T) for T = 0, 1, 2, . . . Tmax-i as defined above, and then extrapolates this monotone 
sequence by a variety of standard sequence extrapolation methods. A similar idea was used to 
furnish estimates of in [3]. 



^We thank John Cardy for this observation. 
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y 

Figure 3: Square lattice with surface site interactions. Axixcy) versus y for T - 1 . . . 15. Inset 
shows the intersection region in finer scale. 



In Section [3] we describe the derivation of the generating functions AT{xc-,y) by the finite- 
lattice method for a range of strip widths T that are needed in this study. For the value of the 

critical step fugacity Xc, we use the exact result Xc = l/\/2 + \/2 for the honeycomb lattice, and 
the best available series estimates in the case of the square and triangular lattices. These are 
Xc{sq) = 0.37905227776 [HI [13], with uncertainty in the last digit, and Xc(tr) = 0.2409175745 
[15] , with similar uncertainty. We performed a sensitivity analysis of our critical surface fugacity 
estimates in order to determine how sensitive they are to uncertainties in our estimates of Xc- 
The estimates of Xc are sufficiently precise that a change in our estimate of Xc by a factor of 10 
times the estimated uncertainty will not change our estimates of the surface fugacity yc in even 
the least significant digit. 

In Section|4]we estimate the critical fugacity by extrapolating yc{T) using a range of standard 
extrapolation algorithms. These are Levin's w-transform, Brezinskii's 6 algorithm, Neville tables, 
Wynn's e algorithm and the Barber-Hamer algorithm. Descriptions of these algorithms, and 
codes for their implementation, can be found in However, we find the most precise estimates 
are given by the Bulirsch-Stoer algorithm [3]. This algorithm requires a parameter w, which 
can be thought of as a correction-to-scaling exponent. For the purpose of the current exercise, 
we have set this parameter to 1, corresponding to an analytic correction, which is appropriate 
for the two-dimensional SAW problem [S] . Our implementation of the algorithm is precisely as 
described by Monroe [TO, and we retained 50 digit precision throughout. 

We used this method to estimate the critical fugacity for all cases of interest for two- 
dimensional SAWs. For the honeycomb lattice, discussed in Subsection 14.11 we have already 
proved [2] that ?/c = 1 + \/2 for the alternate site interaction model, as conjectured by Batchelor 
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and Yung [T] . It is a straightforward consequence of this resuh - the argument is given in Sub- 
section |4?1] below - that for the honeycomb lattice with surface edge interactions (rather than 

site interactions), the critical fugacity is \/l + \/2. For the honeycomb lattice site interaction 
problem where every surface site interacts with the walk, we find the critical fugacity to be 
ijc = 1.46767 where the error in this estimate (and all such estimates given below), is expected 
to be confined to a few parts in the last quoted digit. We know of no other estimate of this 
quantity in the literature. 

In subsection 14.21 we discuss the critical fugacity for site and edge weighted adsorption on 
the square lattice. The only previous estimates for the site weighted case can be found in [12], 
where Monte Carlo methods were used to obtain the estimate Hc = 1.76 ± 0.02. Our estimate, 
t/c = 1.77564 is three orders of magnitude more precise than this. For the edge weighted case, a 
transfer matrix estimate is given in [8j, and is yc = 2.041 ± 0.002. In (7j a Monte Carlo estimate 
of comparable precision is given, yc = 2.038 ± 0.002. Our estimate is yc = 2.040135, again some 
three orders of magnitude more precise. 

For the triangular lattice, discussed in Subsection 14.31 we are unaware of any previous in- 
vestigations of the critical fugacity. We find, in Section 14.31 that ydsite) = 2.144181 and 
j/c(edge) = 2.950026. We repeat that errors in our quoted estimates are expected to be con- 
fined to a few parts in the last quoted digit. 

3 Enumeration of self-avoiding walks 

The algorithms we use to enumerate SAW interacting with a surface on the honeycomb, square 
and triangular lattices builds on the algorithm outlined in our previous paper [3] and detailed 
descriptions can be found in these papers [211121 [H]. Suffice to say that the generating functions 
for a given strip were calculated using transfer matrix (TM) techniques. The most efficient im- 
plementation of the TM algorithm generally involves bisecting the finite lattice with a boundary 
(this is just a line in our case) and moving the boundary in such a way as to build up the lattice 
site by site. If we draw a SAW and then cut it by a line we observe that the partial SAW to 
the left of this line consists of a number of loops connecting two edges in the intersection, and 
at most two unconnected or free edges. The other end of the free edge is an end-point of the 
SAW, hence there are at most two free ends. 

The sum over all contributing graphs is calculated as the boundary is moved through the 
lattice. For each configuration of occupied or empty edges along the intersection we maintain 
a generating function Gs for partial walks with configuration S. In exact enumeration studies 
Gs would be a truncated two-variable polynomial Gs{x,y) where x is conjugate to the number 
of steps and y to the number of surface-contacts (sites or edges) . In a TM update each source 
configuration S (before the boundary is moved) gives rise to a few new target configurations S' 
(after the move of the boundary line) and n = 0, 1 or 2 new edges and m = or 1 new contacts 
are inserted leading to the update Gs'{x,y) = Gs'{x,y) + x"y™Gs{x,y). Here we are primarily 
interested in the case where A(x,y) or B{x,y) are evaluated at the critical point x = Xc- This 
actually makes life easier for us since we can change to a single variable generating function 
Gs{y) and update signatures as Gs'{y) = Gs'{y) + x'2y"'Gs{y)- Here Gs{y) is a polynomial 
in the contact fugacity y with real coefficients truncated at some maximal degree M. The 
calculations were carried out using quadruple (or 128-bit) fioating-point precision (achieved in 
FORTRAN with the REAL(KIND=16) type declaration). 

In our calculations we truncated A{xc,y) at degree M = 1000 and used strips of half-length 
L = M. In Table[T]we have listed estimates for ydd) obtained from strips of width 10 and 9 (the 
crossing between AiQ{xc,y) and Ag{xc,y)) for various values of M and L. Clearly the choice 
M = L = 1000 suffices to estimate yd^) to more than 10 digits accuracy. 

The transfer-matrix algorithm is eminently suited for parallel computations and here we 
used the approach first described in [13] and refer the interested reader to this publication for 
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Table 1: The estimated value of yc(9) from the crossing between Aiq{xc, y) and Ag{xc, y) truncated 
at degree M and using strips of half-length from M up to lOM. 



M 



L^2M L = 5M L = lOM 



100 
250 
500 
1000 



1.832547814756 
1.776250937231 
1.775990340341 
1.775990291271 



1.778376701255 1.778024722094 1.778024722094 
1.775990603337 1.775990594686 1.775990594686 
1.775990291271 1.775990291271 



further detail. The bulk of the calculations for this paper were performed on the cluster of the 
NCI National Facility, which provides a peak computing facility to researchers in Australia. The 
NCI peak facility is a Sun Constellation Cluster with 1492 nodes in Sun X6275 blades, each 
containing two quad-core 2.93GHz Intel Nehalem CPUs with most nodes having 3GB of memory 
per core (24GB per node). It took a total of about 3300 CPU hours to calculate AxixcV) for T 
up to 15 on the square lattice. It is known [14) that the time and memory required to obtain the 
number of walks in a strip of width T grows exponentially as 3"^ for the honeycomb and square 
lattices and as 4^ for the triangular lattice. So, the bulk of the time was spent calculating A15 
and i?i5, which amounted to almost 2300 hours in the square lattice case. In this case we used 
48 processors and the split between actual calculations and communications was roughly 2 to 
1 (with quite a bit a variation from processor to processor). Smaller widths can be done more 
efficiently in that communication needs are lesser and hence not as much time is used for this 
task. 

4 Data analysis 
4.1 Honeycomb lattice 

In [2] we proved that the critical fugacity for the case of interactions with alternate sites on the 
honeycomb lattice is j/c = 1 + ^/2. There are two other cases to consider. The first is the case of 
interactions with every surface site, and the second is the case of interactions with every edge. 
We will deal with the second case first, as it is a straightforward consequence of the proof given 

in [2] that jjc = Vl + in the second case. The proof of this result, in outline, is the following: 
We denote the generating functions A and B, as defined in Section [21 for the alternate site 
case considered in [2], by subscript a (for alternating). We denote the corresponding generating 
functions for the case with edge weighting with the subscript e. Then it is clear by inspection 
that Ac{xc,y) = Aa{xc,y'^), as every time a walk contributing to the A generating function 
passes through n alternating surface sites, whether adjacent or not, it must pass through 2n 
surface edges. 

By the same argument, every time a walk contributing to the B generating function passes 
through n alternating surface sites, whether adjacent or not, it must pass through 2n- 1 surface 
edges. This then gives rise to Bc{xc,y) = ^Bii{xc,y^)- From either of these two equations it 

follows that i/c( alternating) = (yc(edge))^, hence ydedge) = \/l + \/2. 

We now consider the first case, in which every surface site carries a fugacity y. We generated 
data for At{xc, y) for T < 14 as described in Section[3l and found the intersection points where 
Arixcy) = AT+i{xc,y), which defines ydT). These data are tabulated in Tabled Extrapolat- 
ing yc{T) as described above, we estimate 



yc = 1.46767. 



7 



We also find, by an identical method of extrapolation, that A(xc, j/c) = 2.613, which is probably 
exactly 1/ cos(37r/8), as is the case when considering interactions with every alternate site, 
see dH). 

Table 2: The value of yc{T) estimated from the crossing of AT{xc-,y) and AT+i{xc.,y) for the 
honeycomb lattice surface site model. 



T 


yc{T) 


AT{xc,yc{T)) ^ AT+l{Xc,yc{T)) 


1 


1.474342684974343 


2.758023465753132 


2 


1.471231066324457 


2.699581979117133 


3 


1.469859145369675 


2.671309655463187 


4 


1.469144651946551 


2.655387366045945 


5 


1.468728339703417 


2.645467247042683 


6 


1.468465540675101 


2.638829094236329 


7 


1.468289428840316 


2.634145423791235 


8 


1.468122140755486 


2.629489693948282 


9 


1.468008309717543 


2.626054066036805 


10 


1.467956382495343 


2.624432487387554 


11 


1.467915603443970 


2.623117304368586 


12 


1.467883002922926 


2.622033892173660 


13 


1.467856536243392 


2.621129346334020 



4.2 Square lattice 

We next consider data for the square lattice, with every surface site (vertex) carrying a fugacity y. 
We generated data for At{xc, y) for T < 15 as described in Section|3l and found the intersection 
points where AxiXcV) = AT+i{xc,y), which defines ydT). These data are tabulated in Table|31 
Extrapolating ydT) as described above, we estimate 

yc = 1.77564. 

We also find, by an identical method of extrapolation, that A{xc,yc) = 2.678405, which is 
1.024981/ cos(37r/8). In [3] we found, for the non-interacting case (corresponding to y = 1), 
A{xc, 1) = 2.678365 = 1.024966/ cos(37r/8). Thus there appears to be a very weak y dependence. 
(In the normalization of the generating function AT{xc,y) used here, two extra half-steps are 
included, giving an extra factor of the step fugacity Xc, compared to the value that would be 
quoted if contributing walks started and ended on the surface. This explains the difference 
between the values quoted in Table |3] and the ordinates in Figure [31 ) 

Table |4] shows the corresponding data for the edge-weighted case. Extrapolating yc{T) as 
described above, we estimate 

y^ = 2.040135. 

We also find that A{xc,yc) = 2.678405, which is 1.024981/ cos(37r/8). In 3 we found, for the 
non-interacting case (corresponding to y = 1), A{xc, 1) = 2.6783 = 1.0249/ cos(37r/8). This is too 
imprecise to see any evidence of y dependence. 

4.3 Triangular lattice 

We next consider data for the triangular lattice, with every surface site (vertex) carrying a 
fugacity y. We generated data for Axixcy) for T < 11 as described in Section [3l and found the 
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Table 3: The value of ydT) estimated from the crossing of At^XcV) and At+i^XcU) for the 
square lattice surface site model. 



T 


VciT) 


AriXr.VriT)) = AT+^iXr,VriT)) 


1 


1.781782909906119 


2.748677355944862 


2 


1.778386591113354 


2.715115253913871 


3 


1.777378005442640 


2.704018907440273 


4 


1 77fi8^n4n7nQ3'^fi4 




5 


1.776527700942633 


2.693512738663579 


6 


1.776316359764735 


2.690608915840792 


7 


1.776170974231462 


2.688500944397294 


8 


1.776066934443028 


2.686918847615982 


9 


1.775990033953699 


2.685698355993929 


10 


1.775931645420429 


2.684735010917280 


11 


1.775886299456907 


2.683959815456866 


12 


1.775850398954429 


2.683325675630414 


13 


1.775821502307431 


2.682799521958416 


14 


1.775797906369155 


2.682357553489197 



Table 4: The value of yc{T) estimated from the crossing of At{xc, y) and At+i{xc, y) for the square 
lattice surface edge model. 



T 


yc{T) 


AT{xc,yc{T)) = AT^i{xc,yc{T)) 


1 


2.023317607727152 


2.519464246890523 


2 


2.031649211433080 


2.585125356952430 


3 


2.035085448834840 


2.616332757155513 


4 


2.036771224259312 


2.633293109539552 


5 


2.037723730407517 


2.643677266387231 


6 


2.038317002192238 


2.650588857893349 


7 


2.038712823877066 


2.655469267857106 


8 


2.038990695898482 


2.659069610531442 


9 


2.039193569770578 


2.661816780067225 


10 


2.039346383471084 


2.663969985883853 


11 


2.039464457297598 


2.665695001241074 


12 


2.039557641399558 


2.667102372510593 


13 


2.039632511102958 


2.668268404182947 


14 


2.039693596208206 


2.669247312794744 
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intersection points where At{xc, y) = At+i{xc, y), which defines yc{T). These data are tabulated 
in Table [3l Extrapolating yc{T) as described above, we estimate 

j/c = 2.144181. 

We also find, by an identical method of extrapolation, that A{xc,yc) = 4.97002, which is 
1.901944/ cos(37r/8). In 3 we found, for the non-interacting case (corresponding to y = 1), 
A{xc, 1) = 4.970111 = 1.901979/ cos(37r/8). Thus there again appears to be a very weak y depen- 
dence. 

Table 5: The value of yc{T) estimated from the crossing of At^Xcv) and AT+i{xc-,y) for 
triangular lattice surface site model. 



T 


VciT) 


^T^Xc^yciT)) ^ AT+iixc,yc{T)) 


1 


2.169017975620833 


5.299883162257977 


2 


2.152124186067447 


5.089804987842667 


3 


2.147952081330057 


5.033100087535114 


4 


2.146325209334416 


5.009022287728647 


5 


2.145537862947824 


4.996485228732837 


6 


2.145102964455591 


4.989109337635192 


7 


2.144840361941141 


4.984402909686655 


8 


2.144671215263562 


4.981219362650799 


9 


2.144556764080381 


4.978968525942606 


10 


2.144476246964690 


4.977320728801566 



Table m shows the corresponding data for the edge weighted case. Extrapolating yc{T) as 
described above, we estimate 

yc = 2.950026. 

We also find that A{xc,yc) = 4.9696, which is 1.90178/ cos(37r/8). In [3; we found, for the non- 
interacting case (corresponding to y = 1), A{xc,l) = 4.970111 = 1.901979/ cos(37r/8). Again, 
there is evidence of weak y dependence. 

5 Conclusion 

We have estimated the critical fugacity for surface adsorption for two-dimensional SAW on all 
regular lattices for both the case of site and edge interactions. Many of these estimates are new. 
Those that are not are several orders of magnitude more precise than pre-existing estimates. 
Uniquely for the case of the honeycomb lattice with edge interactions, we give the exact value 
of the critical fugacity, and also prove it. Our results are summarised in Table [T] 
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Table 6: The value of VciT) estimated from the crossing of Aj'{xc-,y) and AT+i^x^y) for the 
triangular lattice surface edge model. 



1 


UcK-L ) 


^T{Xc,yc{-l- )) - AT+l{Xc,yc{-l- )) 


-1 
1 











A RA^ 99QR1 0097^4'^ 


3 


2.942788011875285 


4.873934294210283 


4 


2.944814166604381 


4.895179517868169 


5 


2.946090146548846 


4.909648090189844 


6 


2.946944189466541 


4.919989731979732 


7 


2.947544335340955 


4.927679988442194 


8 


2.947982663246637 


4.933582932189477 


9 


2.948312910101248 


4.938231892866670 


10 


2.948568146735367 


4.941971526310544 



Table 7: Estimated critical fugacity yc for surface adsorption. 



Lattice 


Site weighting 


Edge weighting 


Honeycomb 


1.46767 


\/l + V2 


Square 


1.77564 


2.040135 


Triangular 


2.144181 


2.950026 
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